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The low-lying spectrum of the Dirac operator is predicted to be universal, within three classes, 
depending on symmetry properties specified according to random matrix theory. The three universal 
classes are the orthogonal, unitary and symplectic ensemble. Lattice gauge theory with staggered 
fermions has verified two of the cases so far, unitary and symplectic, with staggered fermions in the 
fundamental representation of SU(3) and SU(2). We verify the missing case here, namely orthogonal, 
with staggered fermions in the adjoint representation of SU(A^c), Nc = 2,3. 
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Random matrix theory (RMT) has been successful 
in predicting spectral properties of QCD-like theories 
in the so-called microscopic scaling regime, defined by 
1/Aqcd << L << I/ttItt, with L the linear extent of 
the system |ll01- Assuming spontaneous chiral symme- 
try breaking in the underlying theory, this is the regime 
dominated by the soft pions associated with the chiral 
symmetry breaking S . Up to a scale, given by the chiral 
condensate (at infinite volume) S = (V'V')j the distri- 
bution of the low lying eigenvalues, the so-called micro- 
scopic spectral density 



Ps(z) = lim — pi I 



(1) 



with p{\) the usual (macroscopic) spectral density and, 
in particular, the rescaled distribution of the lowest eigen- 
value Pmin(z) with z = FEAmin are universal, depen- 
dent only on the symmetry properties, the number of 
dynamical quark flavors and the number of exact zero 
modes, i.e., the topological sector, but not the potential 
in RMT B. Recently, these properties have been de- 
rived directly from the effective, finite-volume partition 
functions of QCD of Leutwyler and Smilga, without the 
detour through RMT g. 

The lattice fermion action should have a chiral symme- 
try for the predictions of random matrix theory to apply. 
Until recently, staggered fermions were the only fermions 
regularized on the lattice that retained a chiral symme- 
tryQ The RMT predictions have been nicely verified for 
staggered fermions in the fundamental representations of 
SU(2) II and SU(3) § gauge group in quenched QCD, 



*A new lattice regularization of massless fermions with good 
chiral properties, and even an index theorem, has recently 
been developed H. The RMT predictions for these overlap 
fermions have been verified for examples in all three ensem- 
bles, including the classification into different topological sec- 
tors, in ItI. 



and for SU(2) also with dynamical fermions ||l^. These 
represent two out of the three different cases predicted 
by RMT, chiral symplectic and chiral unitary. Staggered 
fermions are not always in the same universality class 
as continuum fermions in the context of random matrix 
theory. Staggered and continuum fermions belong to the 
unitary ensemble when the fermions are in the fundamen- 
tal representation of SU(A^c) doe Nc > 3. But staggered 
fermions in the fundamental representation of SU(2) be- 
long to the symplectic ensemble since their entire spec- 
trum is two-fold degenerate jO] in contrast to continuum 
fermions which belong to the orthogonal ensemble [||. 
Another example where the staggered fermions and con- 
tinuum fermions are not in the same universality class 
is when the fermions are in the adjoint representation of 
SU(A'c). Here the situation is just reversed: continuum 
fermions are in the symplectic ensemble with a two- 
fold degeneracy of the entire spectrum^] while staggered 
fermions are in the orthogonal ensemble - they are real, 
since adjoint gauge fields are real and since the "Dirac 
matrices" for staggered fermions are just phases, ±1. We 
therefore use this example with staggered fermions in the 
adjoint representation of SU(iVc) to test the missing case. 
In this report, we consider staggered fermions in 
the adjoint representation of SU(2) and SU(3) in the 
quenched approximation. Since staggered fermions (at 
finite lattice spacing) do not obey an index theorem and 
thus have no exact zero modes [Of, the RMT predic- 
tions for ^ = are expected to apply. This has been 
found in the studies with staggered fermions in the fun- 
damental representation p-pll, and our results will sup- 
port it for staggered fermions in the adjoint representa- 
tion. The rescaled distribution of the lowest eigenvalue, 
z = yUAinin, should thus be |0 



^To compare with RMT predictions only one of each degen- 
erate pair of eigenvalues is kept and the adjoint fermion is 
thereby considered as a Majorana fermion. 



TABLE I. The chiral condensate, E, from fits of the distri- 
bution of the lowest eigenvalue to the RMT predictions. The 
last column gives the confidence level of the fit. 



Group 


P 


L 


iV 


s 


Q 


SU(2) 
SU(2) 
SU(2) 
SU(2) 
SU(2) 
SU(2) 
SU(2) 


1.8 
2.0 
2.0 
2.2 
2.2 
2.4 
2.4 


4 
4 
6 
4 
6 
6 
8 


5000 
5000 
2000 
5000 
2000 
2000 
1000 


1.796(18) 
1.676(17) 
1.722(30) 
1.489(16) 
1.529(28) 
1.212(21) 
1.264(30) 


0.921 
0.640 
0.496 
0.840 
0.344 
0.852 
0.874 


SU(3) 
SU(3) 


5.1 
5.1 


4 
6 


2500 
1500 


4.724(76) 
4.671(91) 


0.786 
0.035 



spectrum of low lying eigenvalues of the staggered Dirac 
operator in the adjoint representation very well for both 
SU(2) and SU(3) in our quenched simulations. 
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This distribution is quite distinct from all other cases: it 
starts at a finite, non-zero value at z = 0. Given the chi- 
ral condensate, E, (Q) is a parameter free prediction. If S 
is not otherwise known, it can be determined from a one- 
parameter fit of the distribution of the lowest eigenvalue 
to %. 

We have computed the low lying eigenvalues of the 
staggered Dirac operator in the adjoint representation 
with the Ritz functional method (l4|, applied to D^ D — 
—D^, for several SU(2) gauge field ensembles and for two 
SU(3) ensembles. The distribution of the lowest eigen- 
value, approximated by a histogram with jackknife er- 
rors, was fit to the predicted form, (0), with S, the in- 
finite volume value of (V'V') ^-s the only free parameter. 
The results of the fit, together with the number of config- 
urations considered, gauge coupling and lattice size, are 
given in Table |. In all cases we obtained good fits to the 
predicted form. For most gauge couplings we considered 
two different lattice sizes. The values for E obtained from 
the two different lattice sizes agree well, within statisti- 
cal errors. Some distributions, together with the fitted 
analytical predictions, are shown in Figure H. 

The consistency of the extracted value for the chiral 
condensate, E, from the two different lattice sizes makes 
it evident that we could have used the value obtained 
from one lattice size to get a parameter free description of 
the results from the other lattice size. Given the scale E 
the microscopic spectral density ps, eq. (P, is predicted 



by results of RMTJ 
tems as in Figure [l 



The comparison for the same sys- 
is shown in Figure 0. The agreement 
is quite nice, extending over the entire region covered by 
the eigenvalues that we determined. Note that for the 
chiral orthogonal ensemble the oscillations in the RMT 
prediction for ps, except the first one coming from the 
lowest eigenvalue, are very small. Obviously, we would 
need much more statistics to resolve additional "wiggles" 
in our data. 

In conclusion, adjoint staggered fermions are argued 
to belong to the chiral orthogonal ensemble of RMT. 
We found that the RMT predictions indeed describe the 



[1] 

[2] 
[3] 

[4] 
[5] 



[6] 

[7] 

[8] 



[9] 

[10] 

[11] 
[12] 



[13] 
[14] 



[15] 



E. Shuryak and J.J.M. Verbaarschot, Nucl. Phys. A560 
(1993) 306. 

J. Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531. 
H. Leutwyler and A. Smilga, Phys. Rev. D46 (1992) 
5607. 

G. Akemann, P.H. Damgaard, U. Magnea and S. Nishi- 
gaki, Nucl. Phys. B487 (1997) 721; M.K. Sener and 
J.J.M. Verbaarschot, Phys. Rev. Lett. 81 (1998) 248. 
P.H. Damgaard, Phys. Lett. B424 (998) 322;; G. Ake- 
mann and P.H. Damgaard, Nucl. Phys. B528 (1998) 
411; and Phys. Lett. B432 (1998) 39 0; J.C. Osborn 
D. Toublan and J.J.M. Verbaarschot, tiep-th/980611C| ; 
P.H. DamKaa rd_, J.C. Osborn , D. Toublan and J.J.M. 
Verbaarschot 



hep-th/9811212 



H. Neuberger, Phys. Lett. B417 (1998) 141. 

R.G. Edwards, U.M. Heller, J. Kiskis and R. Narayanan, 



hep-th/9902117 



M.E. Berbenni-Bitsch, S. Meyer, A. Schafer, J.J.M. Ver- 
baarschot and T. Wettig, Phys. Rev. Lett. 80 (1998) 
1146; J.-Z. Ma, T. Guhr and T. Wettig, Eur. Phys. J. 
A2 (1998) 87. 

P.H. Damgaard, U.M. Heller, A. Krasnitz, Phys. Lett. 
B445 (1999) 366; M. Goc keler, H. Hehl, P .E.L. Rakow, 

A. Schafer and T. Wettig, [hep-lat/9811018 . 

M.E. Berbenni-Bitsch, S. Meyer and T. Wettig, Phys. 

Rev. D58 (1998) 071502. 

M.A. Halasz and J.J.M. Verbaarschot, Phys. Rev. Lett. 

74 (1995) 3920. 

J.C. Vink, Phys. Lett. B210 (1988) 211; J.B. Kogut, 

J.F. Lagae and D.K. Sinclair, Nucl. Phys. B (Proc. 

SuppL) 63 (1998) 433; A.L. Kaehler, Nucl. Phys. B 

(Proc. SuppL) 63 (1998) 823. 

P.J. Forrester, Nucl. Phys. B402 (1993) 709. 

B. Bunk, K. Jansen, M. Liischer and H. Simma, DESY- 
Report (September 1994); T. Kalkreuter and H. Simma, 
Comput. Phys. Commun. 93 (1996) 33. 

J. Verbaarschot, Nucl. Phys. B426 (1994) 559. 



0.6 - 



0.4 



0.2 



0.0 




- 0.6 - 



0.4 



0.0 



0.2 



0.0 




0.6 



0.4 



0.2 



0.0 




0.6 



0.4 



0.2 



0.0 




FIG. 1. The distribution of the rescaled lowest eigenvalue, Pminiz), together with the RMT prediction using the best 
estimate for the chiral condensate. 
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FIG. 2. The microscopic spectral density, ps, together with the RMT prediction |l5[ using the best estimate for the chiral 
condensate. 



